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Abstract 


The power law and reflection emission have been observed in the X-ray spectra of both black hole X-ray binaries 
(BHXRBs) and active galactic nuclei (АСМ), indicating a common physical origin of the X-ray emission from 
these two types of sources. The relevant parameters describing the shape of both components and the potential 
correlation between these parameters can provide important clues on the geometric and physical properties of the 
disk and the corona in these sources. In this work, we present a positive correlation between the photon index Г 
and the reflection strength R for the low-mass BHXRBs in the hard state by modeling NuSTAR data, which is 
qualitatively consistent with the previous studies. We compare our results with the predictions from different 
theoretical disk-corona models. We show that the R — Г correlation found in this work seems to favor the moving 
corona model proposed by Beloborodov. Our results indicate that the coronal geometry varies significantly among 
BHXRBs. We further compare our results with that of AGNs. We find that the reflection strength R is smaller than 
unity in the hard state of BHXRBs, while it can be as large as ~5 in AGNS, which implies that the variations of the 
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disk-coronal geometry of AGNs are more vigorous than that of the BHXRBs in the hard state. 
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1. Introduction 


Black hole X-ray binaries (BHXRBs) and active galactic 
nuclei (AGNs) are powered by the accretion of matter onto Ше 
central stellar-mass black holes and supermassive black holes, 
respectively. The accretion matter forms a disk around the 
black hole and liberates the gravitational energy, which is 
radiated in the form of electromagnetic waves (mostly thermal 
emission if the mass accretion rate is not very low, Shakura & 
Sunyaev 1973). The thermal emission from the disk is mainly 
in the optical and UV bands for AGNs and in the X-ray band 
for BHXRBs. In addition, a non-thermal component, which 
cannot be produced in thermal accretion disk accreting at 
20.01 Eddington scaled rate (Meyer et al. 2000b; Maccarone 
2003; Dexter et al. 2021), is also observed in both AGNs and 
BHXRBs. It is widely believed that the soft photons can Бе 
inverse Compton scattered in a region of hot electrons around 
the compact object, namely the corona, in BHXRBs and AGNs, 
which produces the non-thermal primary hard X-ray emission 
(e.g., Galeev et al. 1979; Haardt & Maraschi 1991, 1993; 
Cao 2009; You et al. 2012). However, the geometry of the 
corona is still under debate both observationally and theoretically 
(e.g., Dauser et al. 2013; Liu & Qiao 2022; Yang et al. 2022). It 
has been proposed that the corona could be a radiatively 
inefficient hot accretion flow at the innermost region of the disk 
(Esin et al. 1997; Poutanen et al. 2018), a slab of hot plasma over 


the disk (Petrucci et al. 2001; Wilkins & Fabian 2012), 
a compact hot plasma above the black hole (Miller et al. 2013; 
Ghosh & Laha 2021), or even the base of jet (Markoff et al. 
2003; Shidatsu et al. 2011; Wang et al. 2021). The geometry of 
the corona could also change (e.g., Kara et al. 2019; Méndez 
et al. 2022) during outburst. 

Except for the thermal emission and primary X-ray emission, 
a reflection component has also been observed in both stellar 
and supermassive black hole accretion systems, which is 
produced by the illumination of the cold material by the 
primary X-ray emission, including the fluorescent emission 
lines and a hump in ~20-40 keV energy range (Fabian et al. 
1989; Young et al. 1999; Ross & Fabian 2005; García et al. 
2020). The reflected emission produced at a region close to the 
black hole may also be affected by relativistic effects 
(Laor 1991; Fabian et al. 2000; Reynolds & Nowak 2003; 
Dauser et al. 2012). 

Generally speaking, the primary X-ray emission can be well 
described by a power law with a photon index Г. A cutoff at the 
high energy band has also been reported in some sources with 
high-quality broad-band X-ray data (Fabian et al. 2015; Molina 
et al. 2019; Dong et al. 2022). The reflection strength R 
(К = 0/27, Q is the solid angle covered by the reflector viewed 
from the corona) represents the fraction of the primary X-ray 
photons intercepted by the disk. Despite the great achievements 
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in understanding the accretion processes from the study of 
BHXRBs, some fundamental questions have still remained 
unanswered. For instance, the evolution of the inner region of 
the accretion disk (e.g., Plant et al. 2014; Stiele & Kong 2017; 
Connors et al. 2022), the geometric and physical properties of 
the corona (e.g.,Malzac et al. 2001; Wilkins et al. 2014; You 
et al. 2021; Barua et al. 2022; Méndez et al. 2022), or the 
physical process triggering the accretion state transition (e.g., 
Meyer et al. 2000a; Dong et al. 2022; Liu et al. 2022). These 
questions can be investigated by the study of the evolution of 
specific parameters, e.g., the inner radius of the accretion disk, 
electron temperature, and optical depth of the corona, or the 
correlations between these parameters, e.g., Г and R (Zdziarski 
et al. 1999; Malzac et al. 2001; Molina et al. 2009; Zappacosta 
et al. 2018; Ezhikode et al. 2020), Г and Eddington scaled 
luminosity (Shemmer et al. 2006; Wu « Gu 2008; Gu & 
Cao 2009; Brightman et al. 2013; Yan et al. 2020), inner radius 
and Eddington scaled luminosity (Allured et al. 2013; García 
et al. 2015; Chainakun et al. 2021). 

In this work, we investigate the correlation between the 
photon index Г and reflection strength R in the hard state of a 
sample of stellar-mass black hole systems by modeling their 
Nuclear Spectroscopic Telescope Array (NuSTAR, Harrison 
et al. 2013) energy spectra and compare the correlation in 
BHXRBs with that in AGNs (Panagiotou & Walter 2019; Kang 
et al. 2020; Hinkle & Mushotzky 2021). The correlation 
between I' and R was first found to be positive for the BHXRB 
GX 339-4 by fitting spectra using Ginga observations (Ueda 
et al. 1994). Zdziarski et al. (1999) presented a similar 
correlation using a larger Ginga sample, including Seyfert 
AGNs, BHXRBs, and weakly magnetized neutron stars. The 
R —T correlation was also investigated using RXTE/PCA for 
the BHXRBs GX 339-4 and Cyg X-1 (Revnivtsev et al. 2001; 
Gilfanov et al. 1999). A positive R — Г correlation has only 
been reported for X-ray binaries in the hard state, for instance, 
in Cyg X-1, GX 339-4, GS1354-64, Nova Muscae and GRS 
1915+105 (Gilfanov et al. 2000; Rau & Greiner 2003; 
Zdziarski et al. 2003). 

NuSTAR covers the broad bandpass of 3-79 keV and has 
unprecedented sensitivity above 10 keV, making NuSTAR ideal 
for studies of the disk reflection and the primary emission (Diaz 
et al. 2020; Draghis et al. 2020; Connors et al. 2021). Therefore, 
the launch of NuSTAR provided a new opportunity to study 
К — Г correlation. Using NuSTAR data, ће А — Г correlation 
has been recently confirmed for Seyfert galaxies (Ezhikode et al. 

2020; Panagiotou & Walter 2020) and radio galaxies (Kang 
et al. 2020). Based on the disk-corona model, in which the hot 
electrons in the corona are cooled by the soft photons from the 
cold disk (Haardt & Maraschi 1993), this correlation might be 
caused by the perpendicular movement of the corona with 
respect to the black hole (Beloborodov 1999). If the corona is 
assumed to be moving away from the disk, both reflection and 
irradiation toward the disk are reduced while the reduction of 


Dong, Liu, & Cao 


the soft photons into the corona leads to the harder primary 
emission (a smaller photon index), and vice versa. In Qiao & 
Liu (2017), they proposed that this observed correlation found 
in Seyfert galaxies can be roughly explained within the 
framework of the condensation of the hot corona onto the cold 
accretion disk with the movement of the inner radius of the disk. 

In this paper, we present the sample of BHXRBs and AGNs 
in Section 2. Section 3 gives the spectral analysis. We present 
the best-fitting results and discussions in Section 4. Summaries 
are shown in Section 5. 


2. Sample 


We studied the NuSTAR energy spectra of the BHXRBs 
presented in Yan et al. (2020). Their sample includes 17 
LMXRBs with 165 NuSTAR observations. The data were 
downloaded from HEASARC.? They performed data reduction 
using NUPIPELINE, a task of the NUSTARDAS included in 
the HEAsoft 6.25, with calibration files of version 20181030. 
Further, the source and background spectra were generated 
using NUPRODUCTS. The spectra were extracted from a 
circular region within 90" at the source position. The 
background spectra were extracted from an annulus region 
with inner and outer radii of 180" and 200", respectively. 

In order to achieve the goals of this work, we selected data as 
follows. (1) We excluded observations in which the spectra are 
dominated by the backgrounds (background subtracted photon 
counts below zero). (2) Here, we define the hardness ratio (HR) 
as the 10—78 keV count rate divided by the 3-10 keV count 
rate. To select observations that are in the hard state, we also 
excluded observations with HR<0.3.* (3) The observations in 
which a reflection component is not required in the X-ray 
spectral fitting are also excluded. (4) We further excluded 
observations that cannot be well fitted (х2 > 2) with our 
model (constantxTBabs (diskbb+Gausstpexrav) ог 
constantxTBabs (Gausstpexrav), see the detail in 
Section 3). The final sample consists of 69 observations in 
total. We present the details of the observations in Table 2. 

In order to investigate the possible specific features 
appropriate for black holes of different scales, we will compare 
BHXRBs’ results with different types of AGNs’ in Section 4. 
AGN is classified as Seyfertls (Seyfert 1-1.9) and Seyfert2s 
(Seyfert 2) based on the emergence of the broad emission lines 
in the optical spectrum, or the radio-quiet and radio-loud AGNs 
based on the ratio of the radio to optical luminosity. The 
Seyfert sample is from Hinkle & Mushotzky (2021), including 
26 Seyfertls and seven Seyfert2s, respectively. The radio-quiet 
AGN sample including 28 sources and the radio-loud AGN 
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^ Seein Table 2, there may be some observations (the photon index is ~2.22 
in obs id 80401312002 and ~2.1 in obs id 90202045002 for GRS 1915+105) 
are not in the hard state, but most of the observations (more than 9796 with 
photon index <2.1) are most likely in the hard state (Remillard & 
McClintock 2006). 
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Table 1 

The Column Density of Black Hole Systems 
Source Мн (10? ст?) References 
GRS 1716-249 0.6 Bharali et al. (2019) 
GRS 1739-278 2 Miller et al. (2015) 
GRS 19154-105 6 Feroci et al. (1999) 
GS 1354-64 2 Dunn et al. (2010) 
GX 339-4 0.3 Hynes et al. (2004) 
H 1743-322 22 McClintock et al. (2009) 
IGR J17091-3624 1.2 Rodriguez et al. (2011) 
MAXI J1348-630 0.64 Zhang et al. (2020) 
MAXI J1535-571 5.5 Dong et al. (2022) 
MAXI J1813-095 1.3 Kennea et al. (2018) 
MAXI J1820+070 0.15 Uttley et al. (2018) 
Swift J1753.5-0127 0.15 Steiner et al. (2016) 


sample including 45 sources are from Kang et al. (2020) and 
Panagiotou & Walter (2019), respectively. We refer to these 
papers for detailed information on the observations of 
these AGNs. 


3. Data Analysis 


Spectral fitting is performed using XSPEC 12.11.1 
(Arnaud 1996) with ж statistics. The errors are given with 
9096 confidence level. For each observation, the spectra of 
NuSTAR/FPMA and NuSTAR/FPMB in 3-78 keV are jointly 
fitted, accounting for a cross-normalization with a multi- 
plicative constant We use model TBabs to calculate the 
Galactic absorption. The Verner et al. (1996) photoelectric 
cross sections along with the Wilms et al. (2000) set of solar 
abundances are adopted. Because only fitting the energy band 
above 3keV is not effective to well constrain the column 
density Ny, we fix the Мн to the values in Table 1. 

To model the primary X-ray emission and the reflection 
component, we use the model pexrav, which consists of a 
power-law with an e-folded and a reflection component 
(Magdziarz & Zdziarski 1995). In this model, we assume that 
the abundance of heavy elements, including Fe, of the 
reflecting medium is unity. For the inclination angle i of the 
disk, there are only five sources with constraints on i (60? + 5? 
in GRS 1915+105 (Reid et al. 2014), <79° in GS 1354-64 
(Casares et al. 2009), 37?—78? in GX 339-4 (Heida et al. 2017), 
63° + 3? in MAXI J1820+070 (Atri et al. 2020), and 240? in 
Swift J1753.5-0127 (Neustroev et al. 2014)) in the sample of 
BHXRBs studied in this work. The і cannot either be well 
constrained for the majority of sources. Therefore, for 
simplicity, the inclination angle г is fixed at the model's 
default value, i.e., cosi = 0.45 (i 63°). The photon index Г, 
reflection parameter R, high energy cut Еси, and normalization 
are free parameters. We also include a Gaussian line with 
centroid energy fixed at 6.4 keV to model the Fe ko line. Even 
if the centroid energy is allowed to vary between 6.4 and 
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6.7 keV, the best-fitting results are not further improved. The 
width and the normalization of the line are free. However, the 
width of the line could not be constrained by the data when it is 
not significant enough to be detected, and then, it is fixed at 
0.05 keV. A possible disk component (diskbb) will be 
included if it can significantly improve the fitting. We note that 
the popular reflection model, e.g., ге1хі11 (Dauser et al. 
2014; García et al. 2014), can self-consistently model the 
continuum reflection component and the Fe line (Dong et al. 
2020a, 2020b; García et al. 2022b). However, the pexrav 
model is used in this work so that we can directly compare our 
results with previous studies of NuSTAR AGN samples, in 
which the pexrav model is widely adopted (Panagiotou & 
Walter 2019; Kang et al. 2020; Hinkle & Mushotzky 2021). 

Deviations between FPMA and FPMB of NuSTAR data at 
the lowest energies in the observations for MAXI J1348-630 
(MJD 58515-58516, and MJD 58655-58672) and MAXI J1820 
4-070 (MJD 58192-58397) are found. We let the temperature 
kT and the normalization parameter of the diskbb untie 
between the two detectors for these observations. Note that the 
deviation at the low energy (Е < ~ 6.0 keV) between ЕРМА 
and FPMB in the observations of GX 339-4 on MJD 58060 is 
even worse. Therefore, for this observation, we fit the spectrum 
within 6.0—78.0 keV. The adopted model provides statistically 
accepted fits to our observations (х2 < 2) We list the 
temperature of the disk (KT), the photon index (Г), the 
reflection parameter (R), and the width of the Gaussian line 
(Fe ka с) of the best-fiting in Table 2. The spectra fits are 
insensitive to the Fe abundance, except for GX 339-4 and 
MAXI 71820--070. In order to check the effect of Fe 
abundance on ОХ 339-4 and MAXI J1820+070, we refit 
spectra of the two sources with Fe abundance allowed to be 
free. For GX 339-4, the value of the Fe abundance has a minor 
effect on the parameters К and Г. Therefore, we only report the 
results with Fe abundance fixed at unity for GX 339-4. For 
MAXI J1820-070, a stronger negative correlation between R 
and Г is shown when the Fe abundance is free, which is beyond 
the scope of this paper and will not be further discussed. 


4. Results and Discussions 


This work presents the first systematic study on NuSTAR 
X-ray spectra of BHXRBs to investigate ће А — Г correlation. 
In this sample, we only study the spectra obtained when the 
sources are in the hard state. The finally selected spectra consist 
of 69 observations in total. We modeled the continuum 
reflection component and the Fe line with pexrav and 
Gauss. The photon index Г and reflection strength R are well 
constrained. The values of Г in most observations (more than 
9790) are smaller than 2.1, which is consistent with the 
properties of hard state observation (Remillard & McClintock 
2006). The Г approximately equal to 2.2 and 2.1 in 
observations of GRS 1915+105 (Obs 19 is 80401312002 and 


Table 2 
NuSTAR Observation Details and Spectral Fitting Results 
Source Obs id Detector Exp. Counts HR kT Г R Ее Ка o x и 
(s) (871) (keV) (keV) 
GRS1716-249 80201034006 FPMA 6046.16 113.70 + 0.14 0.47 1.68706? 0.417005 1.431026 2364.31/2309 
FPMB 6134.09 108.00 + 0.13 0.46 
GRS1716-249 80201034007 FPMA 42,461.98 129.20 + 0.06 0.48 0.641020 1.671003 0.381005 1.52 + 0.17 3717.26/3460 
FPMB 43,511.86 119.30 + 0.05 0.47 
GRS1716-249 90202055002 FPMA 17,895.89 124.80 + 0.08 0.40 0.677019 1.791004 0.39+0:06 1.640:15 2960.76/2824 
ЕРМВ 18,258.12 115.00 + 0.08 0.39 
GRS1716-249 90202055004 FPMA 15,796.07 121.10 + 0.09 0.38 0.741008 1.82+005 0.34 + 0.08 151501 2672.78/2649 
ЕРМВ 16,099.26 112.00 + 0.08 0.36 
GRS1739- 278 80002018002 FPMA 29,713.64 110.60 + 0.06 0.35 1.55 + 0.01 0.32 + 0.03 0.99 + 0.04 3032.81 /2474 
FPMB 30,408.71 103.50 + 0.06 0.34 
GRS1915+105 10002004001 FPMA 14,695.91 193.90 + 0.12 0.32 1.92 + 0.01 0.46 + 0.04 1.27 + 0.05 2892.73 /2395 
FPMB 15,151.08 181.10 + 0.11 0.32 
GRS1915+105 80401312002 FPMA 26,166.49 82.62 + 0.06 0.32 2.22 + 0.01 0.92 + 0.05 1.69 + 0.04 2912.76/2484 
ЕРМВ 26,511.82 76.05 + 0.05 0.32 
GRS1915+105 90202045002 FPMA 13,902.54 168.50 + 0.11 0.31 2.13 + 0.01 0.714007 1.59 + 0.06 2522.07 /2374 
FPMB 14,205.90 157.40 + 0.11 0.30 
GRS1915+105 30502008002 FPMA 25,402.94 10.22 + 0.02 0.91 0.99 + 0.03 3.451026 0.45 + 0.01 2695.16/1869 
ЕРМВ 25,434.70 9.44 + 0.02 0.88 
GS1354-64 90101006002 FPMA 24,021.91 6.90 + 0.02 0.60 0.25*016 1.481002 0.28 + 0.06 0.050 1783.01/1656 
ЕРМВ 24,039.58 6.33 + 0.02 0.58 
GS1354-64 90101006004 FPMA 28,799.04 47.93 + 0.04 0.52 0.521010 1.541004 0.641012 1.411049 2778.47/2626 
ЕРМВ 28,921.35 44.38 + 0.04 0.50 
GS1354-64 90101006006 FPMA 34,968.15 54.60 + 0.04 0.47 0.521010 1.591002 0.431004 1.197013 2892.15/2720 
FPMB 35,032.79 51.51 + 0.04 0.46 
GX339- 4 80001013002 FPMA 35,093.22 10.36 + 0.02 0.53 1.59 + 0.01 0.31 + 0.04 0.38016 2071.64/2081 
FPMB 35,101.11 10.10 + 0.02 0.52 
GX339-4 80001013004 FPMA 20,796.25 17.29 + 0.03 0.53 1.59 + 0.01 0.36 + 0.04 0.36010 2224.51 /2072 
FPMB 20,836.18 16.30 + 0.03 0.52 
GX339-4 80001013006 FPMA 23,154.34 27.56 + 0.03 0.52 1.61 + 0.01 0.39 + 0.03 0.481009 2767.26/2425 
ЕРМВ 23,246.65 25.91 + 0.03 0.51 
GX339-4 80001013008 FPMA 43,193.17 36.17 + 0.03 0.51 1.64+001 0.45 + 0.02 0.551007 3875.44/2960 
FPMB 43,442.76 34.28 + 0.03 0.50 
GX339-4 80001013010 FPMA 97,091.69 6.13 + 0.01 0.50 1.60 + 0.01 0.23 + 0.03 0.38005 2428.24/2361 
FPMB 97,785.98 5.84 + 0.01 0.49 
GX339-4 80001015001 FPMA 10,376.54 103.20 + 0.10 0.38 1.94 + 0.01 0.76 + 0.07 1.407007 2417.03 /2324 
FPMB 10,549.13 95.83 + 0.10 0.37 
GX339-4 80102011002 FPMA 18,464.23 18.65 + 0.03 0.39 1.92 + 0.01 0.621007 1.03 + 0.11 2215.44/1781 
ЕРМВ 18,194.09 17.40 + 0.03 0.38 
GX339-4 80102011004 FPMA 16,849.29 16.33 + 0.03 0.45 1.77 + 0.01 0.46 + 0.06 0.96015 2011.26/1781 
FPMB 17,384.07 15.21 + 0.03 0.44 
GX339-4 80102011006 FPMA 16,462.32 13.80 + 0.03 0.47 [71509 0.38007 0.81703 1721.02/1693 
FPMB 16,579.54 12.94 + 0.03 0.45 
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Table 2 
(Continued) 
Source Obs id Detector Exp. Counts HR kT Г R Ее Као x V 
(s) (s 1) (keV) (keV) 
GX339- 4 80102011008 FPMA 19,933.55 10.80 + 0.02 0.47 s 1.68 + 0.01 0.31 + 0.05 0.63016 1798.84/1681 
ЕРМВ 20,197.23 10.25 + 0.02 0.47 
GX339-4 80102011010 FPMA 37,366.39 7.66 + 0.01 0.49 Е 1.66 + 0.01 0.29 + 0.04 0.41795 2004.22/1864 
FPMB 37,185.80 7.22 + 0.01 0.47 
GX339-4 80102011012 FPMA 41,335.80 4.13 + 0.01 0.48 = 1.65 + 0.02 0.221006 0.632016 1569.84/1557 
ЕРМВ 41,326.41 3.84 + 0.01 0.46 
GX339-4 80302304002 FPMA 21,137.86 2.41 + 0.01 0.48 "m 1.64 + 0.03 0.20010 0.29028 1059.93/1017 
FPMB 21,155.48 2.23 + 0.01 0.46 
GX339-4 80302304004 FPMA 17,100.77 21.91 + 0.04 0.51 a 1.62 + 0.01 0.43700 0.747033 2297.29/2061 
FPMB 16,957.45 20.48 + 0.04 0.50 
GX339-4 80302304005 FPMA 18,707.55 19.64 + 0.03 0.51 Е 1.51 + 0.02 0.167005 0.432018 1911.08/1882 
ЕРМВ 18,673.39 18.30 + 0.03 0.49 
GX339-4 80302304007 FPMA 28,964.39 3.72 + 0.01 0.48 а: 1.65 + 0.02 0.25 + 0.07 0.411030 1358.57/1361 
ЕРМВ 28,890.22 3.52 + 0.01 0.47 
GX339-4 80502325008 FPMA 22,397.06 10.32 + 0.02 0.43 = 1.76 + 0.02 0.457008 1.021020 1734.00/1652 
ЕРМВ 22,249.86 9.56 + 0.02 0.42 
GX339- 4 90401369004 FPMA 3621.47 13.03 + 0.06 0.51 = 1.59 + 0.03 0.3670:18 0.321030 987.79/1002 
ЕРМВ 3616.32 11.98 + 0.06 0.50 
H1743-322 80001044002 FPMA 50,257.48 26.10 + 0.02 0.57 Е 1.55 + 0.01 0.347003 1.01704? 3027.51/2874 
FPMB 50,392.11 25.31 + 0.02 0.55 
H1743-322 80001044004 FPMA 61,284.19 34.70 + 0.02 0.54 а 1.58 + 0.01 0.397005 0.924033 3263.14/3052 
FPMB 61,345.92 33.77 + 0.02 0.53 
H1743-322 80001044006 FPMA 25,669.39 30.89 + 0.04 0.55 па 1.57 + 0.01 0.35004 0.947051 2367.63/2484 
FPMB 25,724.12 29.06 + 0.03 0.53 
H1743-322 80002040002 FPMA 27,419.12 17.92 + 0.03 0.56 gs 1.53 + 0.01 0.22 + 0.03 0.497055 2137.19/2250 
FPMB 27,604.23 16.45 + 0.02 0.54 
H1743-322 80202012002 FPMA 65,862.57 38.45 + 0.02 0.45 па 1.75 + 0.01 0.42700 1,827502, 3141.24/2944 
FPMB 66,276.43 36.22 + 0.02 0.43 
H1743-322 80202012004 FPMA 64,763.73 36.83 + 0.02 0.44 gi 1.78 + 0.01 0.427005 1.821021 3063.71/2895 
ЕРМВ 65,033.95 34.46 + 0.02 0.43 
H1743- 322 80202012006 FPMA 65,702.66 20.08 + 0.02 0.54 а 1.57 + 0.01 0.261005 1.031036 3130.68 /2849 
ЕРМВ 65,610.73 18.94 + 0.02 0.53 
H1743-322 90401335002 FPMA 38,447.51 32.20 + 0.03 0.51 s 1.60 + 0.01 0.35 0.03 122702) 2881.08/2700 
FPMB 38,392.55 30.31 + 0.03 0.51 
IGRJ17091-3624 80001041002 FPMA 43,293.43 15.72 + 0.02 0.52 e 1.691006 0.721023 2.912028 2394.47/2365 
FPMB 42,953.51 14.63 + 0.02 0.51 
IGRJ17091-3624 80202014002 FPMA 20,237.91 20.52 + 0.03 0.50 Е 1.63 + 0.01 0.41 + 0.05 0.721025 1969.75/1973 
ЕРМВ 20,276.72 19.08 + 0.03 0.48 
IGRJ17091-3624 80202014004 FPMA 20,698.10 22.66 + 0.03 0.48 I: 1.62 + 0.01 0.32 + 0.04 0.541022 1930.61/1985 
ЕРМВ 20513.13 21.07 + 0.03 0.46 
MAXIJ1348-630 80402315002 FPMA 3038.04 1043.00 + 0.59 0.43 0.50 + 0.07 1.79 + 0.01 0.67 *003 1.62 + 0.11 3183.00/2994 
FPMB 3208.38 904.50 + 0.53 0.42 
MAXIJ1348-630 80402315004 FPMA 736.34 1048.00 + 1.20 0.42 0.5070 1.807002 0.747013 1.827021 2220.32/2187 
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Table 2 
(Continued) 
Source Obs id Detector Exp. Counts HR kT Г R Fe Као x и 
(s) (871) (keV) (keV) 
FPMB 778.54 910.50 + 1.08 0.42 
MAXIJ1348- 630 80502304002 FPMA 13,783.85 186.80 + 0.12 0.45 0.61 + 0.05 1.71 + 0.01 0.39 + 0.03 1.344048 3336.50/3046 
FPMB 14,172.58 171.60 + 0.11 0.44 
MAXIJ1348-630 80502304004 FPMA 15,367.24 161.90 + 0.10 0.46 0.73 + 0.06 1.68 + 0.01 0.36 + 0.03 121700 3232.84/3054 
FPMB 15,773.25 149.00 + 0.10 0.45 
MAXIJ1348-630 80502304006 FPMA 17,181.65 107.50 + 0.08 0.47 0.56 + 0.10 1.68 + 0.01 0.38 + 0.04 1.381025 3008.45 /2923 
FPMB 17,460.58 99.65 + 0.08 0.46 
MAXIJ1535-571 90301013002 FPMA 10,143.05 607.00 + 0.25 0.41 0.44 + 0.03 1.84+0.01 0.79 + 0.04 1.70 + 0.05 3740.87 /3072 
FPMB 10,671.15 546.50 + 0.23 0.40 
MAXIJ1813-095 80402303002 FPMA 20,544.76 16.10 + 0.03 0.49 ves 1.68 + 0.01 0.33 + 0.04 0.05(f) 1974.03/1913 
FPMB 20,538.99 15.05 + 0.03 0.48 
МАХ 1813-095 80402303004 FPMA 20,431.99 13.37 + 0.03 0.49 ot 1.67+0.01 0.30 + 0.05 0.391020 1778.12/1803 
ЕРМВ 20,344.93 12.60 + 0.03 0.48 
МАХІЈІ813-095 80402303006 ЕРМА 23,236.11 14.45 + 0.03 0.46 I 1.75 + 0.01 0.377008 0.731026 1914.72/1854 
ЕРМВ 23,096.79 13.55 + 0.02 0.45 
SwiftJ1753.5-0127 30001148002 FPMA 40438.46 12.28 + 0.02 0.42 e 1.73 + 0.01 0.11 + 0.03 0.05(f) 2101.11/1987 
FPMB 40,480.61 11.62 + 0.02 0.41 
MAXIJ1820--070 90401309002 FPMA 11,768.40 177.30 + 0.12 0.55 1.41 + 0.07 1.44+001 0.24 + 0.02 0.26 + 0.04 3612.24/3207 
ЕРМВ 11,979.89 168.10 + 0.12 0.54 
MAXIJ1820+070 90401309004 FPMA 2760.54 728.70 + 0.52 0.49 0.951009 1.712001 0.971008 2.041015 3398.69/2985 
ЕРМВ 2876.64 657.30 + 0.48 0.48 
MAXIJ1820--070 90401309006 FPMA 4539.49 754.10 + 0.41 0.49 1.177001 1.74 + 0.01 1.347016 2.661019 4075.04/3235 
ЕРМВ 4761.95 673.90 + 0.38 0.48 
MAXIJ1820--070 90401309008 FPMA 3045.78 769.60 + 0.51 0.48 0.89+007 1.75 + 0.01 1.02 + 0.08 2.011014 3490.94/3021 
FPMB 3213.80 682.20 + 0.46 0.47 
MAXIJ1820+070 90401309010 FPMA 2660.13 775.30 + 0.55 0.48 0.89 *007 1.741001 1.04 + 0.09 2.043043 3491.23/2962 
FPMB 2801.20 689.10 + 0.50 0.47 
MAXIJ18204-070 90401309012 FPMA 12,333.18 688.70 + 0.24 0.47 0.87 + 0.02 1.69*009 0.70 + 0.02 1.72 + 0.05 5602.26/3572 
FPMB 12,963.19 617.10 + 0.22 0.46 
MAXIJ1820+070 90401309013 FPMA 1834.23 666.80 + 0.61 0.45 0.807007 1.731001 0.691008 1.771017 2821.99 /2630 
ЕРМВ 1933.90 591.40 + 0.56 0.45 
MAXIJ18204-070 90401309014 FPMA 9208.05 671.40 + 0.27 0.46 0.80 + 0.02 1.72 + 0.01 0.68 + 0.03 1.74 + 0.06 4684.80/3446 
ЕРМВ 9707.75 596.90 + 0.25 0.45 
MAXIJ1820--070 90401309016 FPMA 13,791.26 567.80 + 0.20 0.43 0.83 + 0.02 1.737000 0.55 + 0.02 1.74 + 0.04 4628.45 /3513 
FPMB 14,536.45 509.20 + 0.19 0.42 
MAXIJ18204-070 90401309018 FPMA 2733.29 500.40 + 0.43 0.41 0.832000 1.751001 0.502005 1.611009 2745.80/2616 
ЕРМВ 2879.46 451.90 + 0.40 0.41 
MAXIJ1820+070 90401309019 FPMA 9443.58 499.00 + 0.23 0.42 0.82 + 0.02 1.74 + 0.01 0.51 + 0.03 1.717005 3927.23/3292 
FPMB 9937.79 450.80 + 0.21 0.41 
MAXIJ18204-070 90401324002 FPMA 6570.42 482.00 4- 0.27 0.41 0.83 + 0.03 1.74 + 0.01 0.46 + 0.03 1.681007 3597.19/3104 
ЕРМВ 6931.00 432.80 + 0.25 0.41 
MAXIJ18204-070 90401309021 FPMA 21,867.84 299.20 + 0.12 0.43 0.81 + 0.02 1.64 + 0.01 0.31 + 0.02 1.50 + 0.04 3973.26 /3460 
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Table 2 
(Continued) 
Source Obs id Detector Exp. Counts HR kT D R Ее ka o x и 
(s) (571) (keV) (keV) 
FPMB 22,852.97 270.40 + 0.11 0.42 
MAXIJ1820+070 90401309033 FPMA 24,894.75 123.10 + 0.07 0.31 0.81 + 0.01 1.89+001 0.35 + 0.03 1.74 + 0.03 3072.78 /2858 
FPMB 25,345.13 115.40 + 0.07 0.31 
MAXIJ1820+070 90401309035 FPMA 18,571.50 52.83 + 0.05 0.44 1.73 0.01 0.307005 0.753029 3042.14/2526 
FPMB 18,642.35 50.28 + 0.05 0.43 
MAXIJ18204-070 90401309037 FPMA 37,682.99 14.05 + 0.02 0.42 1.74 + 0.01 0.17 + 0.03 0.537313 2060.92 /2047 
FPMB 37,660.36 13.15 + 0.02 0.41 
MAXIJ1820+070 90501311002 FPMA 28,664.18 13.68 + 0.02 0.42 1.75 + 0.01 0.24 + 0.04 0.317002 1926.63/1891 
FPMB 28,588.09 12.80 + 0.02 0.42 
MAXIJ18204-070 90501337002 FPMA 44,577.56 14.24 + 0.02 0.43 1.71 € 0.01 0.20 + 0.03 0.28 + 0.05 2327.92/2184 
ЕРМВ 44,429.89 13.27 + 0.02 0.43 
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Figure 1. The best-fitting parameters Г апа R. The gray lines are calculated based on the model of moving corona (Beloborodov 1999). js represents the geometry of 


the corona, which corresponds to и, = 0, 0.1, 0.2, 0.3, 0.4, and 0.5. 


Obs id is 90202045002, respectively), which may be asso- 
ciated with the model, will not affect the conclusion in this 
paper. For MAXI J1820-F-070, the Г is in the narrow range of 
~1.6-1.8, while Ше R spans a wide range of —0.1—1.4. These 
two parameters in MAXI J18204-070 do not follow the positive 
R — Г correlation reported in previous works. We also note that 
the best-fitting parameter Г is 0.9 + 0.03, and R is 345027 in 
the spectrum of GRS 1915+105 on MJD 58623 with 
x = 2695.16 / 1869. We consider these results as outliers of 
the regular R — Г panel, which will not be further discussed in 
our work. For the remained observations, the values of R 
obtained are roughly in the range of 0.1—1.0, and Г ranges from 
~1.4 to ~2.2. We present the points of R and Г in Figure 1. We 
test the correlation between К and Г using the Spearman 
method. The Spearman's coefficient is equal to 0.58 with a p- 
value of 1.13 x 107°, which indicates а significant positive 
К — Г correlation. 

The А – Г correlation can provide important clues on the 
geometric and physical properties of the disk and the corona. 
Zdziarski et al. (1999) reported an extremely strong positive 
correlation between К and Г in the sample of Seyfert galaxies, 
and the observations of their four X-ray binaries also satisfied 
this correlation. Zdziarski et al. (2003) further evaluated the 
reality of this correlation. In the work of Zdziarski et al. (1999), 
the obtained values of А are smaller than 2 and the values of Г 
are 1.4—2.2. The correlation is considered to provide evidence 
that the seed photons for the primary source come from the 
same medium that is responsible for the reflected emission. 


They found that the model with the inner radius movement 
failed to reproduce the data in the observed range of Г and К, 
but only the data in the mid-range. On the other hand, the 
model proposed by Beloborodov (1999), in which the corona 
above the accretion disk moves toward or away from the black 
hole with relativistic bulk velocity, appears quantitatively to 
reproduce the correlation. This scenario is indeed consistent 
with the X-ray observations of an AGN sample, which provides 
evidence of the outflowing corona in AGNSs (Liu et al. 2014). 
Below, we calculate the correlation for a moving primary 
source using the equations of Beloborodov (1999). In their 
model, the reflection strength R is computed as follows, 


в. 2+ 8/20 - ди) 
(1+ By 


where 8 = и/с is the velocity of the primary source, and 
и = cosi (i is the inclination of disk). We assumed 
cosi — 0.45 as fitting the data using pexrav. The spectral 
index Г is calculated as follows, 


(1) 


Газ 2.33(A — 1), (2) 


where 6— 1/6 and 6=1/10 are for BHXRBs and AGNs, 
respectively. A is the amplification factor which can be 
estimated as follows, 


2ү(1 + BPC + up? 


А = ; 
(1 — Фа — а) — (1 + i)? 6?/4] 


(3) 
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where у = (1 — 82) 1/2. а (az20.1-0.2) is the reflection 
albedo, and ju, is the coronal geometry. The parameter a will be 
assumed to be 0.15 when we calculate. The case u, = 0 
represents a slab geometry of the corona, while и, = 0.5 means 
the corona roughly to be a blob with a radius of order its height. 

We note that the geometry of the corona is still under debate. 
Here, we plot the expected curve for jj, = 0, 0.1, 0.2, 0.3, 0.4, 
and 0.5 as gray lines in Figure 1. It can be seen that all lines fail 
to reproduce the observed range of Г and А. Most of our data 
points fall into the area between the two lines with us = 0 and 
Hs = 0.5 which may indicate that the geometry of the corona 
changes. The case that the coronal geometry evolves has been 
reported in previous studies. The work of Kara et al. (2019) 
found that the reverberation lags between the corona and the 
accretion disk in the BHXRB MAXI J18204-070 shorten while 
the profile of the broad Fe line remained remarkably constant 
and explained this case with a reduction in the spatial extent of 
the corona. In the microquasar GRS 1915-105, it is also found 
that the initially extended X-ray corona can become more 
compact and turned into the jet (Karpouzas et al. 2021; García 
et al. 2022a; Méndez et al. 2022) with the evidence that (1) the 
phase lags of the type-C QPO made a transition from hard to 
soft, (2) the two separate correlations are presented between the 
iron-line flux and the total flux, (3) and the QPO frequency and 
the ratio flares evolved in a systematic way. Note that in the 
model proposed by Beloborodov (1999), R — 1 is given for a 
static corona and independent of the coronal geometry. The 
ranges of 8 used to estimate the theoretical lines are larger than 
О for all assuming values of и, indicating that the corona is 
moving away from the black hole. When the corona moves 
away from the black hole, the beaming effect reduces the 
photons irradiating the disk and subsequently the reflection 
emission. Moreover, the cooling effect on the corona temper- 
ature by seed photons reduces, and then a harder spectrum is 
presented. 

In order to investigate the possible specific features 
appropriate for black holes of different scales, we compare 
BHXRBs’ results with different types of AGNs’. In the work of 
Hinkle & Mushotzky (2021), they used XMM-Newton and 
NuSTAR data to analyze 33 AGNs, including 26 Seyfertls and 
seven Seyfert2s. The authors presented the main fundamental 
parameters, such as Ше Г, К, folding energy, black hole mass, 
and Eddington ratio for Seyfert AGNs, and looked for 
correlations between them. Seyfert2s tend to have lower Г 
(harder spectra) than Seyfertls, and the distributions of R for 
Seyfertls and 2s are consistent with each other. They 
recovered the positive R— Г correlation with p-value = 
2.6 x 107° combining Seyfertls and Seyfert2s. The work of 
Kang et al. (2020) measured Г, К, folding energy, and Fe Ka 
EW for 28 radio-loud AGNs, and compared their results with 
45 radio-quiet AGNs’ (Panagiotou & Walter 2019), in both of 
which NuSTAR observations were analyzed. They found that 
the radio-loud AGNs have harder spectra and weaker reflection 
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emission compared with radio-quiet ones. Kang et al. (2020) 
reported a positive R — Г correlation with a very low p-value 
(«10 ^) combining radio-quiet and radio-loud, which con- 
tributed to the moving corona model. 

We compare the best-fitting R and Г for the hard state of 
BHXRBs and the sample of Seyferts, and for the hard state of 
BHXRBs and the sample of radio-loud and radio-quiet in 
Figures 2 and 3, respectively. In both figures, the distributions 
of R and Г are given in the lower right and upper left panels. 
The significant feature is that the reflection strength of 
BHXRBs' hard state is smaller than 1, but it can increase 
largely to 5 in the AGN sample. When we see the distributions 
of R and Г in Figure 2, there are no significant differences for 
the three populations, while in Figure 3, the hard state of 
BHXRBs tends to have flatter spectra than radio-quiet and 
radio-loud AGNs, weaker reflection than radio-quiet АСМ, 
but the distribution of R is consistent with radio-loud AGNs. To 
check for the robustness of the differences in the hard state of 
BHXRBs and radio-quiet/loud AGNs, we perform two non- 
parametric statistical tests, 1.е., Mann-Whitney U test (hereafter 
referred to as the U-test). The U-test for two independent 
samples can be calculated without the requirement of normal 
distribution. The null hypothesis is that there is no difference 
between any two compared data sets, and it is rejected when the 
p-value is below 0.05. We find that the distributions of Г and R 
in BHXRBs are clearly different from that in radio-quiet 
AGNs, where p-value is 3.1 x 1077 and 2.1 x 1074, respec- 
tively. But BHXRBs and radio-loud AGNs have similar Г and 
R distributions with p-value being larger than 0.05 and 0.86, 
respectively. The relationships between R and Г are shown in 
the lower left panels in Figures 2 and 3. Except for the 
theoretical lines predicted by Beloborodov (1999) model for 
ВНХКВв, we also calculate the expected lines with и, = 0 and 
Hs = 0.9 for AGNs. It can be seen that the areas between these 
two predicted lines mostly cover the observational points of 
AGN ignoring its classification. It is likely that the BHXRBs 
need a narrower area than AGNSs to reproduce the data of А and 
Г, which may imply that the corona in BHXRBs changes less 
vigorously than that in the AGNs. 

In the low/hard state of the BHXRBs, an advection- 
dominated accretion flow (ADAF) is present in the inner 
region near the BH (Narayan & Yi 1994), extending out to 
connect a thin disk at a truncated radius (Esin et al. 1997). The 
truncated radius Ry of the thin disk (i.e., the ADAF radial size) 
increases with the decrease in the mass accretion rate (e.g., 
Yuan & Narayan 2014). Thus, the reflection strength R may 
decrease with increasing truncated radius R,. As only the 
observations in the hard state of the BHXRBs have been 
included in our present investigation, while the AGN samples 
employed in this work include some luminous AGNs, which 
roughly correspond to the BHXRBs in the high/soft state. The 
inner edge of the thin disk in those luminous AGNs may extend 
to the inner stable circulating orbit (ISCO), and therefore their 
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Figure 2. Distributions of photon index Г and reflection strength R in the upper left and lower right panels, and R vs. Г in the lower left panel. The results of Seyferts 
are from Hinkle & Mushotzky (2021). The orange, magenta, and black are for Seyfertls, Seyfert2s, and BHXRBs, respectively. The gray solid and dotted lines 
correspond to 44, = 0 and us = 0.5 for BHXRBs, and the blue solid and dotted lines correspond to ptis = 0 and и, = 0.9 for AGNs. 


Fraction 
о о 
н N 


Èo 
o 


o 
л 


0.0; 


-3.0 


BHXRB: и; = 0 
BHXRB: и; = 0.5 
AGN: р; = 0 
AGN: us = 0.9 
Radio loud 
Radio quiet 
BHXRB 


1.00 


1.25 


1.50 


2.25 2.50 0.0 0.1 0.2 0.3 0.4 0.5 0.6 


Fraction 


Figure 3. Distributions of photon index Г and reflection strength R in the upper left and lower right panels, and R vs. Г in the lower left panel. The results of radio-loud 
are from Kang et al. (2020), and radio-quiet are from Panagiotou & Walter (2019). The orange, magenta, and black are for radio-loud, radio-quiet, and BHXRBs, 
respectively. The gray solid and dotted lines correspond to ps = 0 and jj, = 0.5 for BHXRBs, and the blue solid and dotted lines correspond to ps = 0 and ps = 0.9 


for AGNs. 
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reflection strength К is higher than that of the BHXRBs. It is 
also worth noting that the outer radius of a thin disk may vary 
from several thousand to >10 Schwarzschild radii (You et al. 
2012; Cao 2016; Cao & Lai 2019; Cao & Zdziarski 2020), 
however, the affection of the outer disk radius to the reflection 
strength К should be much less important compared with that 
of the truncated radius of the thin disk. Thus, we think that it 
contributes little to the difference of R between BHXRBs 
and AGNs. 


5. Summary 


We investigate the spectra of the hard state in black hole 
low-mass X-ray binaries using NuSTAR data. A positive 
correlation between К and Г is tested by the Spearman method, 
with the Spearman's coefficient equal to 0.54 and a p-value 
equal to 4.27 x 107°, except for the source MAXI J1820+070. 
К and Г in MAXI J1820+070 seem to be the outlier in the 
regular R — Г panel, which will be further explored in a future 
paper. The К-Г correlation presented here has been reported in 
GX339-4 and Cyg X-1, which is consistent with the previous 
report for AGNs. Our results support the moving corona model 
(Beloborodov 1999), in which the corona moves perpendicu- 
larly with respect to the black hole. To mostly cover the whole 
observational points, the geometry of the corona needs to be 
changed. When we compare the best-fiting R and Г for the 
BHXRBs and the AGNs, we find that Ше А is smaller than 1 in 
BHXRBs while it can be as large as ~5 in AGNs, which 
implies that the variations of the corona in AGNs are more 
vigorous than that of the BHXRBs in the hard state. We 
conjecture that the luminous AGNS with thin discs extending to 
the ISCOs may also be responsible for the systematically 
higher R in AGN samples. 
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